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Abstract 

A good understanding of electroluminescence is a prerequisite when optimising double-phase 
noble gas detectors for Dark Matter searches and high-pressure xenon TPCs for neutrinoless 
double beta decay detection. 

A simulation toolkit for calculating the emission of light through electron impact on neon, 
argon, krypton and xenon has been developed using the Magboltz and Garfield programs. Calcu- 
lated excitation and electroluminescence efficiencies, electroluminescence yield and associated 
statistical fluctuations are presented as a function of electric field. Good agreement with experi- 
ment and with Monte Carlo simulations has been obtained. 

Keywords: Electroluminescence, Electron drift. Noble gases. Dark matter, Neurinoless 
double-beta decay. 



1. Introduction 

Several experiments in astrophysics and cosmology, such as direct Dark Matter searches lU]- 
[stl and neutrinoless double beta decay [4, 5] , are based on noble gas and/or liquid Time Projec- 
tion Chambers (TPC) which use electroluminescence for primary ionisation signal amplification. 
This technique gives high gains and good energy resolution, and is suitable for experiments with 
low event rates and high background levels. 

Up to now, for calculating the light yield and efficiencies, a three-dimensional Monte Carlo 
program of the electron drift in xenon and xenon-neon mixtures 01 and a one-dimensional 
program for krypton and argon [8,] existed. Although validated, they are not open-source nor 
freely accessible, unlike the toolkit described in this paper Our toolkit calculates the yield and 
efficiencies and provides a versatile and comprehensive parameterisation of the physics processes 
involved in electroluminescence. It uses an up-to-date database of electron-atom and electron- 
molecule cross sections for about 60 gases. 

In this paper, calculations for neon, argon, krypton and xenon in uniform fields are discussed. 
We compare our results with earlier Monte Carlo simulations j^-lst] and with experimental data 
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Figure 1 : Simplified energy level diagram of argon. The energy levels were taken from 111511 and 
the ionisation thresholds from [16]. The stronger dipole-allowed transitions are also shown lll7ll . 
Continuous arrows (purple) correspond to VUV, dash-dotted arrows (red) to visible and dashed 
arrows (brown) to IR photons. 



2. Electroluminescence 

In a gaseous detector, primary electrons are produced along the track of a incoming particle 
interacting in the absorption region. They are driven by an electric field, below the scintillation 
threshold, to a region where, under the influence of stronger fields, they can excite or ionise the 
gas atoms. 



The excited atomic levels are described by the jl coupling H12LI13I1 . In the Racah notation, 
e.g. the lowest four excited levels are ns [3/2]2, ns [3/2]°, ns' [1/2]q and ns' [1/2]°. Fig lT]shows 
the level diagram of argon with the stronger dipole-allowed transitions ][l3]- 

In addition to the transitions shown in Fig. [T] experimental spectra of pure noble gases 
show a continuum in discharges and in proportional scintillation 1 181420(1 due to excimer decay. 
Excimers - electronically excited molecules, - are formed through three-body collisions 
between an excited atom, R*, and two atoms in the ground state, R: 

R*+2R^R"+R. (1) 

Three-body collisions dominate at pressures above a few tens of Torr making the excimer 
formation the main channel of de-population of excited atoms [21]. In this paper we assume 
that other processes, e. g. associative ionisation which affects highly excited states that are not 



frequently produced, contribute negligibly Il22i 12311 . 
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The excimers involved are mainly and ^2+ JUB. They are formed through process 
([U in high vibrational states and can decay to the repulsive ground state, 'E^, emitting a VUV 
photon, 



2R + hvu 



(2) 



or they can collide with ground state atoms losing vibrational energy: 



Rl*+R 



Rl+R. 



(3) 



In the latter case, the resultant excimer in a low vibrational state, ' emits a VUV photon 
with slightly lower energy: 



r;^2R + hv2. 



(4) 



The electronic transitions of excimers follow the Franck-Condon principle 11261 12711 . Tak- 
ing into account that the ground state is repulsive, the principle explains the continuum spectra 
observed experimentally for low pressures (< 100 Torr) which exhibit two peaks, the "first con- 
tinuum" at higher frequencies - and the "second continuum" at lower frequencies. 

At high pressures, typically above 300 Torr, the proportional electroluminescence spectra 
show only the second continuum because process Q dominates over process ^ ifisll . 



3. Simulation toolkit 



The simulation toolkit was developed using the new C-n- version of the microscopic tech- 
nique of Garfield [28], which currently uses Magboltz 8.9.3 ll29ll . 



5.7. Garfield 

Garfield is a program for the detailed simulation of gaseous detectors ^3^. Its Monte Carlo 
microscopic technique tracks electrons in gases at molecular level using procedures and cross- 
sections available in Magboltz. 

In noble gases each collision is classified as elastic, excitation or ionisation. 

Information about each excited atom is available: the (x, y, z) position, the time of production 
and the excitation level. This information is used to determine the electroluminescence signal 
properties. 

The program simulates virtually any electric and magnetic field Bill . 



3.2. Magboltz 

Magboltz was developed to calculate the transport parameters of electrons drifting in gases 
under the influence of electric and magnetic fields [32]. 

For this purpose, the program contains, for 60 gases, electron cross sections for all relevant 
interactions. Separate excitation cross sections are available for 45, 44, 4 and 50 levels/groups of 
neon, argon, krypton and xenon, respectively. 

Fig. |2]presents the total number of collisions of different types produced during electron drift 
over z - 2 cm of gas at 293 K and 1 atm - the conditions of our simulations - as a function of 
the reduced electric field, E/N, i.e. the electric field E divided by the gas number density A^. The 
ns states represent ~90 % of all excitations in the proportional electroluminescence region (see 
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Figure 2: The total number of collisions (elastic, ionisations and of some excitation 
groups/levels) during electron drift over z - 2 cm of gas at 293 K and 1 atm, including those 
experienced by secondary charges when ionisation is possible. 
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Figure 3: Mean and maximum energies, e and Smax respectively, reached by the electrons before 
collisions, as functions of EjN. 



Fig. IT) in accordance with fis"]. The np states represent between 5 - 10 %, depending on the gas 
and on the intensity of the field. 

The most frequently produced excited state in Ar andXe is the metastable level ns [3/2]2 with 
a lifetime of seconds |33ll. closely followed by the radiative «s states, with an intrinsic lifetime 
of nanoseconds lf34ll35[l . Atomic transitions from the latters are promptly re-absorbed by ground 
state atoms through the radiation trapping mechanism. The long effective lifetime of these states 
enhances the importance of excimers. 

Fig. [3] shows the maximum energy of electrons, over 10'^ collisions, as well as their mean 
energy. 



4. Model 

We assume that each excited atom gives rise to the emission of one VUV photon from second 
continuum excimer decay. Its wavelength was generated according to a Gaussian distribution 
with mean 82.2 nm, 128 nm, 147 nm, and 173 nm and a FWHM of 3 nm, 10 nm, 12 nm and 
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14 nm for neon, argon, krypton and xenon, respectively 118-23, 36]. In addition, when charge 
multiplication occurs, each ion also gives rise to the emission of one VUV photon iIstIi . 

The number of primary electrons for each EjN was varied between 3x10"^ and 2 x 10^ to 
keep the error in J (see Eq. O) below 2 %. The starting direction of each primary electron was 
sampled isotropically. The starting energy was distributed according to the energy distribution 
calculated by Magboltz. 



5. Results and discussion 

Fig. |4]shows the excitation efficiency, Qexc, the fraction of energy acquired by the A'e primary 
electrons in the electric field that is spent in the excitation process ia]. 



ezN^E 

where «exc is the number of excitation groups available in Magboltz for the studied gas, n' the 
number of excitations of the /* group produced by the primary electrons, e'^^^ the energy of the 
excitation group, z the distance travelled by the electrons, E the electric field and e the elementary 
charge. 

Fig. |5] shows the electroluminescence efficiency, Qel, the ratio between the energy emitted 
in the form of VUV photons and the energy acquired by the electrons during drift: 

i=l i=l 

Qel = ^TT— (6) 

where s'^^ is the energy of the de-excitation photon of the excited atom in the group through 
excimer decay. This parameter was generated randomly according to the characteristics of the 
second continuum (Section|4]i. 

When the electric field supplies enough energy to electrons for excitation, the cross section 
increases abruptly and electroluminescence thus begins with a fast increase of 2exc and Qel- 
Increasing E/N, a plateau is reached before Qexc and Qel increase due to the electroluminescence 
produced by the additional secondary electrons. 

Qel is always lower than Qexc- This is due to the loss of vibration energy from excimers 
before they emit a VUV photon, to radiative transitions from molecular levels to and 
and to infrared losses. 

2exc and Qel in argon and krypton agree for high E/N with f^, 's']. Below 6 to 8 Td 
(1 Td = 10-'^ Vcm^) our values are 10 % lower The difference increases with decreasing E/N. 
This region is not interesting for operating detectors based on electroluminescence amplification 
because of the large statistical fluctuations and the low light gains. 

The first Townsend coefficient a, obtained dividing the number of ionisations per primary 
electron by the drift distance, is shown in Fig. |6] 

Fig. |7] shows the reduced electroluminescence yield Y/N, defined as the number of photons 
emitted per primary electron and per unit of drift path divided by A^. Fig. |7]also presents results 
of earUer simulation work for xenon and measurements for argon and xenon [jgl [l^l. The 
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overall agreement consolidates the assumption that the main channel of de-population of excited 
atoms is through excimer decay. If other processes contribute they also lead to the emission of a 
VUV photon. The exception in the agreement is argon below 8 Td and is under investigation. 

YjN is linear in EjN at low EjN, even when some ionisation is produced (see Fig. |6]l. Above 
~7 Td, ~15 Td, ~18 Td and ~22 Td for neon, argon, krypton and xenon, secondary charges 
producing electroluminescence, change the linear behaviour of YjN. These thresholds are in 
good agreement with those calculated for neon and xenon in Fig. 5 of ifllil . For lower values 
of EjN, the probability of ionisation is too low for changes in the linear behaviour of YjN to 
be detected. The slopes increase from the lighter to the heavier gas and reflect the decrease 
in the minimum energy required to produce one excitation. Xenon is the gas that gives the 
highest electroluminescence gains in the linear region, followed by krypton, argon and neon. 
Extrapolating the electroluminescence yield, we determined the electroluminescence threshold 
to be 1.51 + 0.04 Td, 4.1 ± 0.1 Td, 2.6 + 0.1 Td and 2.9 + 0.1 Td for neon, argon, krypton and 
xenon, respectively, in good agreement with [ 8, ,1 IJ . Neon, although it has the highest minimum 
excitation energy, has the lowest electroluminescence threshold. This is due to neon having the 
smallest elastic cross section of all gases studied. Thus, the electrons achieve higher energies at 
the same field. 

The relative variance in the number of emitted photons J, is defined as 
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y = ^ (7) 
A^vuv 

where cr^u^j is the variance of the number of photons emitted by one primary electron, A^vuv- 
This parameter, shown in Fig. |8] is useful for determine the energy resolution of a detector 
1 3811 . Xenon is the gas that exhibits the lowest statistical fluctuations, followed by krypton, argon 
and neon. For electric fields above the electroluminescence threshold, as the field increases, J 
decreases until the onset of secondary electron production. At this point, J begins to increase 
because the higher fluctuations in the charge gain start to dominate. To the best of our knowledge, 
up to now, there are not consistent values for J published in literature. 



6. Conclusions 

We used the new C++ version of the microscopic technique of Garfield to access information 
on excited atoms produced during the electron drift through the pure noble gases neon, argon, 
krypton and xenon. The simulation relies on procedures and cross-sections available in Magboltz 
8.9.3. We assumed that every excited atom leads to the emission of one VUV photon. We were 
able to calculate for a uniform electric field: the excitation and electroluminescence efficiencies, 
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the electroluminescence yield and the statistical fluctuations as functions of the reduced electric 
field. 

The results obtained for the excitation and electroluminescence efficiencies and for the exci- 
tation thresholds are in good agreement with earlier Monte Carlo work. We also compared our 
results for the reduced electroluminescence yield with measurements for argon and xenon, and 
good agreement was found. 

We simulated the statistical fluctuations associated with electroluminescence. As the reduced 
electric field increases, the statistical fluctuations decrease until the secondary charge fluctuations 
dominate. We confirm that the statistical fluctuations associated to proportional electrolumines- 
cence are lower than those in both charge avalanche multiplication and primary electron cloud 
formation. 
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Figure 8: Relative variance in the number of emitted photons as a function of EjN for 2 cm of 
drift, 293 K and 1 atm. 
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